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ABSTRACT 

We study the amplification of magnetic fields during the formation of primordial halos. The turbu- 
lence generated by gravitational infall motions during the formation of the first stars and galaxies can 
amplify magnetic fields very efficiently and on short timescales up to dynamically significant values. 
Using the Kazantsev theory, which describes the so-called small-scale dynamo - a magnetohydrody- 
namical process converting kinetic energy from turbulence into magnetic energy - we can then calculate 
the growth rate of the small-scale magnetic field. Our calculations are based on a detailed chemical 
network and we include non-ideal magnetohydrodynamical effects such as ambipolar diffusion and 
Ohmic dissipation. We follow the evolution of the magnetic field up to larger scales until saturation 
occurs on the Jeans scale. Assuming a weak magnetic seed field generated by the Biermann battery 
process, both Burgers and Kolmogorov turbulence lead to saturation within a rather small density 
range. Such fields are likely to become relevant after the formation of a protostellar disk and, thus, 
could influence the formation of the first stars and galaxies in the Universe. 

Subject headings: cosmology: early universe - ISM: abundances - ISM: magnetic fields - turbulence 



1. INTRODUCTION 

Magnetic fields play an important role in star for- 
mation. There is increasing evidence that dynam- 
ically important magneti c fields were also present 
in the early Univ e rse (iBaneriee fc JedamzikT 120041: 
Silk fc Laneeri [2005 ISchleicher et all 120101: [Sur et alJ 



20101: IFederrath et all 1201 Ibl: iTurk et al.ll2012D . If this 
turns out to be true, current models for the formation of 
the first stars and galaxies need to be revisited. 
The theory of primordial star formation has gone 
through a change recently. It was previously assumed 
that t he first stars w e re extremely massive a nd iso- 
lated (jAbeTetaD [IMl iBrom^ELario^J^^ How- 



ever, new high-resolution calculations ( Stacv et al.l 



2010; 



2011 



■i gn-res qlr 

Clark et al.]l20nl: ISmith et al.l IMlt IGreif et al.l 

20121) show that the accretion disk of a collapsing pri- 
mordial halo fragments into multiple stars. The inclusion 
of magnetic fields could change this picture again. The 
magnetic pressure can stabilize an accretion disk and, 
depe nding on the fie ld strength, suppress the fragmenta- 
tion (iMachida et al, 2004; Hcnncb clle fc Tevssier 2007t 
Price fc Batd 120071: iPeters et al.l 120111: iSeifried et'alj 
201 It ). So far, there are only a few stud- 
ies on pri mordial star formation that include mag- 
netic fields (|Tan fc BlackmanI [200l iMaki fc Susall200l 
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). The magnetic field is expected 



to have similar effects to those seen in present-day star 
formation, suc h as the launching of winds and jets 
(jMachida et al. 2006, 2008a). The latter eject gas from 
the accretion disk which could otherwise have collapsed 
onto the star. Thus, the star formati on efficiency is re- 
duced, especially for high-mass stars (|Tan fc BlackmanI 
[200l IMachidal [20l0l) . Strong jets can transport matter 
even out of the star-forming halo, l eading to a ma gneti- 
zation of the intergalactic medium (jXu et al.ll2011l) . 
For developing primordial star formation theory further, 
we need to know structure and strength of the magnetic 
fields. The mechanism we suggest for producing strong 
fields is the small-scale dynamo. This magnetohydro- 
dynamical (MHD) process amplifies weak magnetic seed 
fields exponentially by converting kinetic energy from the 
turbulence i nto m agnetic energy. 

iKazantsevI ()1968l ) developed a theory for describ- 
ing the small-scale or turbulent dynamo. The 
main equation of this theory is the Kazantsev 
equation, the eigenvalue of which is the growth 



ent solutions of this equation ( 


Roeachevskii fc Kleeorin 


1997|: 


Kleeorin fc Roecachevskiil 


20111: iSchekochihin et al. 


2002a 


). In this work we use 1 


he solution proposed by 
<es into account different 


Schober et al.l (|2012D. which tal 



types of turbulence. The growth rates of the magnetic 
field obtained in that work are co mparable to previous 
results for Kolmogorov turbulence (jSubramanianl ITQQTI : 
iBrandenburg fc Subramanianll2005l ). 
We model the physical and chemical processes during the 
collapse of a primordial halo to quantitatively determine 
the properties of the small-scale dynamo. We calculate 
the magnetic Prandtl number and the magnetic Reynolds 
number. The latter is compared to the critical magnetic 
Reynolds number for small-scale dynamo action. Fur- 
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Figure 1. The temperature as a function of tlie number density. 
The different lines indicate different initial conditions. The dashed 
green line corresponds to an initial temperature of 10* K, the dot- 
ted red line to 10^ K and the solid orange line to 10^ K. 

thermore, we calculate the growth rate, which depends 
on the magnetic Prandtl number and the hydrodynamic 
Reynolds number. We assume a weak initial magnetic 
field of 10"^" G o n the viscous scale produced by the 
Biermann battery (IBierm annlllQSCII : IKulsrud et al.lfl997l : 
iXu et all 12008V This allows us to determine the evolu- 
tion of the magnetic field strength during the collapse. 
The structure of our study is as follows. In Section [2] 
we review the properties of primordial gas. We present 
our numerical calculation of the chemistry and thermal 
evolution of the gas and discuss the characteristic mag- 
netohydrodynamical quantities. Furthermore, we discuss 
the origin of turbulence and weak magnetic seed fields in 
primordial halos, which are essential for operation of the 
small-scale dynamo. Section|3]concentrates on the small- 
scale magnetic field evolution. We give the basic equa- 
tions of the Kazantsev theory and the resulting growth 
rates of the magnetic energy. We apply these results for 
the small-scale magnetic field to our model for the col- 
lapse of a primordial halo. In Section |4] we present a 
model for the transport of the magnetic energy to larger 
scales. This allows us to calculate the magnetic energy 
on the Jeans scale of the primordial halo. 

2. PROPERTIES OF THE PRIMORDIAL GAS 

2.1. Chemical and Thermal Evolution 

We determine the chemical and thermal evolution of 
gravitationally collapsing primo r dial ga s using the one- 
zone model of iGlover fc SavinI (120091), togethe r with 
a modification implemented by iSchleicher et al.l (|2OO90 
that relates the collapse time to the equation of state. 
Moreover, we have included additional Li"*" ch emistry 
by using the reaction rates from Bovino et al.l 

and HeH"*" chem i stry a ccording to I Bovino et al , ^ 

IGlover fc SavinI (|200 9^ model the chemistry of the gas 
with a chemical network that includes around 30 different 
atomic and molecular species linked by around 400 dif- 
ferent chemical reactions. In our calculations, we use the 
sa me initial chemical abu ndances as in the default model 
in IGlover fc SavinI ()2009D . The elemental abundances of 
helium, deuterium and lithium relative to hydrogen are 
taken to be 0.083 for helium , 2.6 x 10~^ f or deuterium 
and 4.3 x 10'^° for lithium (|Cvburtll200l . The initial 




Figure 2. The fractional abundances of different chemical species 
as a function of the number density. 



density and temperature of the gas were assumed to be 
no = 1 cm^'^ and To — 1000 K, respectively, but we have 
verified that our results are not sensitive to these values. 
In the one-zone model the mass density p evolves as 



dt tff 



(1) 



where tg = {32Gp) is the free-fall time. Moreover, 

evolution is determined by the energy 



the temperature 

equation, 

de 
dt 



dt 



A 



-heat 7 



(2) 



where e is the specific internal energy, p is the thermal 
pressure and Acooi and Aheat are the total cooling and the 
heating rate per unit mass, respectively. The resulting 
temperature evolution is shown in Figure |T] Besides the 
initial temperature of 10^ K, we also show in this figure 
the results for the initial temperatures 10^ K and 10'* K. 
Our calculations result in roughly the same evolution for 
all the initial temperatures after an increase of the den- 
sity of about one order of magnitude. 
In Figure|2l the fractional abundances of H, He, II2, H"*", 
Li"*" and free electrons vary with increasing density in 
our calculations. The abundance of H is constant at low 
densities, but decreases at densities higher than about 
10^" cm"-^ due to the formation of H2. As there is no 
dust in primordial gas, large quantities of H2 are pro- 
duced only at high densities, via the three-body reac- 
tions: 



H-hH-HH^H2 + H, 
H-hH-HHe^H2 + He, 
H + H + H2^H2 + H2. 



(3a) 
(3b) 
(3c) 



For the magnetic properties of the primordial gas the 
abundances of the charged species are especially im- 
portant. They determine for example the conductiv- 
ity, which is calculated in the next section. At densities 
n < 10^ cm"'^, ionized hydrogen is the main positive ion, 
while at higher densities, Li+ dominates. The sharp drop 
in the 11+ abundance at densities n > 10* cm~'^ results 
from the removal of 11+ from the gas by the reaction 



Small-Scale Dynamo in Primordial Star Formation 



3 



chain () Glover fc Savinl[2009l) . 

H2 + H+-^H 
H++e-^H2 + H 



3 +7, 



(4a) 
(4b) 



IS 



2.2. Characteristic Magnetohydrodynamical Quantities 
Viscosity — It can be shown that the kinematic viscosity 

1/2 



1 



kT 



(5) 



if the molecules are assumed to be rigid spheres 
(jChoudhuril |1998[) . Here, n = p/m is the number den- 
sity, k the Boltzmann constant and T the temperature. 
Furthermore, d — J2i 6*^1 is the mean particle diameter 
and m — fiWi the mean mass, is the relative abun- 
dance of the species i, mi and di are the masses and the 
Van-der-Waals diameters, respectively. The temperature 
as well as the abundances of the individual species are 
functions of the number density. 

Diffusimty — For calculating the magnetic diffusivity 
r/ we need the conductivity of the gas. In a plasma 
the three most important contributions to the conduc- 
tivity of a ne utral species (indicated by index n) are 
(jWardle fc Ngl ll9991. 



C \ ^ 



/3n 



■ 0'^ 



(6a) 



(6b) 



(6c) 



as given in lPinto et al.l (|2008[) . The Hall-parameters /3n 
are defined as 



yns — 



(7) 



Here, and rris are the masses of the neutral and the 
charged particles, and are the abundance fractions 
of the species, and {<7v)^^ is the momentum transfer 
rate coefhcient. We take these c oefficients, wh i ch ar e 
functions of the temperature, from lPinto fc GalU (|2008l ). 
where we use the polarisation approximation for Li+ . 
The two dominant effects that lead to the dissipation of 
magnetic energy are the Ohmic resistivity and ambipo- 
lar diffusion. We can neglect the contribution of the Hall 
effect to the resistivity, as here the force acts perpen- 
dicular to the current and, thus, no energy is dissipated 
into heat. We calculate the distributions of the Ohmic 
resistivity and the ambipolar diffusion by 



VOhi, 



47r(T| 



(8b) 



We focus on the most important neutral species H, He 
and H2 and the charged species H"*", e~ and Li"*". For 
each neutral species we calculate the resistivities Tyohm.n 



and 77AD,n- The magnetic field strength B drops out 
in the Ohmic case. Finally, the total Ohmic magnetic 
diffusivity is 77ohm = J2n ^Ohm,n and the total resistivity 
due to ambipolar diffusion is 77AD ~ l/(Sn^AD n)0 

Reynolds Numbers — The hydrodynamic and magnetic 
Reynolds numbers are defined as 



Re: 



Rmi 



VL 

V 

VL 

7 

7? 



(9a) 
(9b) 



where L is the length of the largest turbulent fluctua- 
tions and V the typical velocity on that scale. Notice, 
that we give these numbers on the forcing scale, i.e. the 
Jeans scale, which means L = £j and V = v^. 
For the calculation of the magnetic Reynolds number we 
use the sum of jyohm and 77ad- The resulting Reynolds 
numbers are shown in Figure [3] as a function of the den- 
sity. The critical magnetic Reynolds numbers (jlSp are 
also indicated for the two extreme types of turbulence. 



Magnetic Prandtl Number- 
netic Prandtl number is 

Pm 



The definition of the mag- 



Rm 



V 

v' 



(10) 



We can calculate this quantity by using the equations ([5]) 
and ([8]). In Figure [3] the density dependency of the mag- 
netic Prandtl number is shown for both Kolmogorov and 
Burgers turbulence. For clarification we point out that 
the rapid decrease of the magnetic Reynolds and Prandtl 
number is caused by the dynamo amplification of the 
magnetic field. In the beginning of the collapse Ohmic 
resistivity is the dominant diffusion process. With in- 
creasing magnetic field ?7ad increases proportional to B^ 
(see equations [6] to [8]) and becomes the main process for 
magnetic diffusion. Since Rm and Pm are both propor- 
tional to I/77AD, in the limit where 77AD ^ VOhm, both 
decrease rapidly with increasing magnetic field strength. 
In addition, we tested the influence of varying the initial 
temperature on the evolution of the Reynolds numbers 
and the magnetic Prandtl number. For the initial tem- 
perature ranging from 10^ K to 10"* K we found only 
small variations in the Reynolds numbers and the mag- 
netic Prandtl number, as illustrated in Figure [31 

2.3. Turbulence 

Turbulence is an omnipresent phenomenon in astro- 
physics. Turbulence is observed, for example, in the 
convection zone of stars, in accretion discs, and in the 
interstellar medium. 

For the first star-forming halos considered here, 
we as sume that t u rbule n ce is drive n by virialisa- 
tion (IWise fc Abell [20071 : IGreif et al.l [2008h or by 
accretion of gas into the center of the halos 
(!kI csscn & Hcnncbellc'MSIlIme green fc BurkertJiMol : 
Fcdcrrath et al. 2011a). The presence of turbu- 
lence affects star formation strongly, as the tur- 
bule nt pressure works against the collapse to a 
star (| Vazquez- Semadeni et al.l[l998HMac Low fc KlessenI 

^ From private communication with Daniele Galli. 
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Figure 3. The hydrodynamic and magnetic Reynolds numbers, 
Re and Rm, as well as the magnetic Prandtl numbers, Pm, on 
the current Jeans scale. The numbers are presented as a func- 
tion of the hydrogen nuclei number density n. The solid lines 
represent an initial temperature of 10^ K, the dashed lines 10^ K 
and the dotted lines 10^ K. Moreover, the horizontal lines indi- 
cate the critical magnetic Reynolds number for Kolmogorov and 
Burgers turbulence (^JmJ^it = 107 and Rm^^^^. = 2718) as derived 
in lSchober et ali (2012). The rapid decrease of Rm and Pm from 
the very high starting values is caused by the exponential dynamo 
amplification of the magnetic field. We show the results for Kol- 
mogorov turbulence in the upper plot and the results for Burgers 
turbulence in the lower plot. 



200l iKrumholz fc McKed [200l iMcKee fc Ostrikeii 
2007h . Moreover. iFederrath et all pOllaD show in a 
Fourier analysis that the turbulence is effectively driven 
on the Jeans scale. 

There are different types of turbulence. In this paper 
we concentrate on the two extreme cases, Kolmogorov 
turbulence and highly compressible Burgers turbulence. 
The different types are described in the inertial range by 
the relation between the length scale £ and the velocity 
V on that scale, 

voce'^. (11) 

The exp onent ranges from 1/3 for incompressible tur- 
bulence (lKolmogorovlll94lD to 1/2 for highly compress- 
ible turbulence ( Burgerslll948[ ) . In real astrophysical ob- 
jects we expect the tu rbulenc e in dex 'd to li e between 
these extreme cases (jKritsuk et al.i 120071 : iSchmidt et al.i 



[2008Hl^dCTrath et aLllMoh . 

2.4. Magnetic Seed Fields 

There are different theories that describe the ori- 
gin of weak primordial magnetic fields. The first seed 
fields could already have bee n produced during infla- 
tion. [Turner fc Widrowl (|1988D find that a magnetic field 
Bq « 10~^^-10~^ nG on a scale of 1 Mpc can be pro- 
duced wh en the con f ormal invariance is broken. 
Following iSigl et all ()1997[) . there is also a possibility to 
create a magnetic field during first-order phase transi- 
tions in the very early Universe. They predict a field 
strength Bq « 1G~^° nG from the electroweak phase 
transition and Bq w 10^^^ nG from the QCD phase tran- 
sition on a scale of 10 Mpc. 

Another popular way to generate magnetic fields are 
battery mechanisms. The so-called Biermann battery 
uses the fact that electrons and ions have very different 
masses. If there is a pressure gradient in the plasma the 
particles get accelerated. Due to their smaller masses 
the electrons are more strongly accelerated than the 
ions. This leads to charge separation and an electric 
field is generated. If the electron density Ue is con- 
stant in space the electric field is static, however, if it 
varies in space electric currents are generated, which 
give rise to a corresponding magnetic field. Note, that 
for the Biermann term to be non-zero the cross prod- 
uct of the electron pressure gradient and the electron 
density gradient needs t o be non-zero (|Biermannl 119501 : 
IKulsrud fc Zweibe]|[200l . 

3. MAGNETIC FIELD AMPLIFICATION ON THE VISCOUS 

SCALE 

In this section we analyse the evolution of the small- 
scale magnetic field. We outline the Kazantsev theory, 
which gives the growth rates of the magnetic field on 
the viscous scale. Together with the amplification due to 
gravitational compression and dissipation processes we 
can calculate the resulting small-scale magnetic field evo- 
lution. 

3.1. Small-Scale Dynamo Growth 

Kazantsev Theory — An arbitrary magnetic field B can, 
in general, be separated into a mean component Bq and 
a fluctuating component (5B with 



B = Bo + SB. 



The induction equation 
dB 
'dt 



V X (v X B - r/V X B) 



(12) 



(13) 



describes the time evolution of this field, where v is 
the velocity and rj the magnetic diffusivity. Substitu- 
tion of ([T^ into the induction equation leads to two 
equations: an equation for the larg e-scale field evo- 
lution and the Kazantsev equatio n ([Kazantsev! 119681 : 
[Brandenburg fc Subramanianll200"5l) , which describes the 
small-scale evolution of the field. 

The derivation of the Kazantsev equation is based on the 
assumption that the fiuctuations of the magnetic field 
as well as the fiuctuations of the velocity field are ho- 
mogeneous and isotropic even if the mean fields are not 
isotropic. Furthermore, the fluctuations are assumed to 
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be a Gaussian random field with zero mean and the ve- 
locity fluctuations are assumed to be delta-correlated in 
time. For the simplicity helicity of the magnetic field is 
neglected. With these assumptions the Kazantsev equa- 
tion is 



d r 



U{r)'iP{r) = -F?A(r). 



(14) 



The eigenfunctions of this equation are related to the 
longitudinal correlation function of the magnetic fluctu- 
ations MLir,t) by Ml = l/{r'^y^K^)^{r)e'^^K We caU 
F the growth rate of the small-scale magnetic field. The 
function Kdiff is the magnetic diffusion coefhcient, which 
contains besides the magnetic diffusivity rj also a scale- 
dependent turbulent diffusivity. U is called the "poten- 
tial" of the Kazantsev equation. Both K^iff and U depend 
only on the correlation function of the turbulent velocity 
field and the magnetic diffusivity. The correlation func- 
tion of the turbulent velocity field in turn depends on the 
different types of turbulence. 

With a mode l for the turbulent correlation function, 
iSchober et al.l ()2012[ ) solved the Kazantsev equation (|14|) 
with the WKB-approximation. They found that the 
critical magnetic Reynolds number for dynamo action 
Rrricrit increases with compressibility. The values that 
Rm needs to exceed are 



ao7, 

^2718, 



(15a) 
(15b) 



for Kolm ogorov and Burgers tu rbulence, respectively. 
Moreover, ISchober et al.l |2012|) found different growth 
rates of the magnetic field for different turbulence mod- 
els, with 

(163-304^y 

60 L ^ ' 

in the limit of infinite magnetic Prandtl numbers. Here 
V is the typical velocity on the largest scale of the tur- 
bulent eddies of size L and Re is the hydrodynamical 
Reynolds number. 

In this paper we analyse the two extreme types of turbu- 
lence, Kolmogorov with -d — 1/3 and Burgers turbulence 
with 1? = 1/2. We find in the limit of large magnetic 
Prandtl numbers 



37V 
36 



-Re 



1/2 



F^ = i?e^/^ 

60 L 



(17a) 
(17b) 



For the typical velocity of the largest fluctuations we 
use the sound speed V = V] ^ y/jkT/m, a s the Mach 
number in a primordial halo is roughly one (|Greif et al.l 
|2008[ ). Here 7 is the adiabatic index. We take L to 
be the Jeans length, as this is the effect ive driving scale 
for turbulence in a collapsin g system (jSchleicher et al.l 
I2010t iFederrath eraIll2011aD . Hence we set L ^ £j = 

\J^kT j {Gm^n), where G is the gravitational constant. 
We compare the growth rate of the small-scale dynamo 
F on the viscous scale £iy = l]Re~^^^^~^^^ to the inverse 
free-fall time l/tg = [37r/(32G'TOn)]-i/2. The result is 
shown in Figure |4l In our model the magnetic field on 
the fastest growing scale increases one to three orders of 
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Figure 4. Different characteristic rates on the viscous scale as a 
function of the density. We compare the growth rate of the mag- 
netic field (solid orange line) to the inverse free fall time (dashed- 
dotted gray line) and the dissipation rates (Ohmic dissipation: 
dashed green line, ambipolar diffusion: dotted red line). We show 
the results for Kolmogorov turbulence in the upper plot and the 
results for Burgers turbulence in the lower plot. 

magnitude faster than the halo collapses. Note that the 
dynamo growth is exponential in time. 

Gravitational Compression — The gravitational compres- 
sion due to the collapse of the halo provides additional 
amplification of the magnetic field. As long as the con- 
dition of flux freezing is fulfilled, the magnetic field B 
increases with density like 



B cx 



(18) 



for spherically symmetric collapse. Before the dynamo 
saturates, the amplification by gravitational compression 
is minor compared to the dynamo growth. 

Dissipation — Part of the magnetic energy is converted 
into heat by dissipation processes. The dissipation term 
in the induction equation (|13l) is TyV^B. We consider 
Ohmic dissipation and ambipolar diffusion and approx- 
imate this by rjB/i'^ and dB/dt by BY ohm and BTab 
respectively. We get 



Ohm ' 



Tad ' 



VOhn 



(19a) 
(19b) 



Tohm and Fad are the rates of magnetic energy dissipa- 
tion by Ohmic resistivity and ambipolar diffusion. 
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3.2. Critical Magnetic Reynolds Number 

The dependency of the magnetic Reynolds number Rm 
on the number density is shown in Figure [3] for the two 
extreme types of turbulence. We also indicate the crit- 
ical magnetic Reynolds number for small-scale dynamo 
action Rm„it- One can see that the magnetic Reynolds 
number is larger than Rrucnt at the onset of the col- 
lapse. For densities above roughly 4 cm~^ Rm becomes 
smaller than the critical value in the case of Burgers 
turbulence. For Kolmogorov turbulence Rm becomes 
also smaller than Rmcmt for high densities, which are 
not shown in Figure [3l However, as we see below, at 
this point the dynamo is already saturated on the small 
scale as well as on the large scale. Thus, in the den- 
sity regimes where the small-scale dynamo operates, the 
condition Rm > Rmciit is always fulfilled. 

3.3. Resulting Small-Scale Magnetic Field 

In principle, the magnetic energy density, Eb — 
i?^/(87r), evolves as 



d^B 
dt 



4 dn 
3n dt 



Ohm 



(20) 



where we assume spherically symmetric collapse. 
By solving equation I^U^ numerically we find the evo- 
lution of the magnetic energy density on small scales. 
In Figure [7] we show the resulting growth of the mag- 
netic field strength. As an initial field strength Bq we 
use 10"^*^ G on the viscous scale, which is a conser- 
vativ e value for a field generated b y a Biermann bat- 
tery ()Biermannlll950l:lXu et al.ll2008D . The field strength 
grows extremely rapidly as the density increases. How- 
ever, we cannot trust the whole evolution of the mag- 
netic field exactly as shown in Figure [T] When the 
field has become strong enough, the magnetic Prandtl 
number becomes unity or less (see Figure ^ . Then the 
WKB-approximation breaks down and equations (ITT]) are 
no longer valid. Complementary studies have shown, 
however, that the small-scale dynamo still operates for 
Pm < 1 (Boldvrcv fc Cattaneo"2004': fSchekochihin et al.l 
\2ml J007.; Evinld I2011D ~ although the growth rate 
may decrease by a f actor of a few. We note that 
iBoldvrev fc Cattaneol (|2004f ) find in their studies that 
the critical magnetic Reynolds number increases with de- 
creasing magnetic Prandtl number. Furthermore, we see 
in Figure |4] that the ambipolar diffusion rate becomes 
higher than the growth rate of the magnetic field. In 
this regime, equation (fTH]) is no longer a solution of ([T4|) . 
We expect that the field grows at the rate (|16l) almost 
until saturation, but then decreases and the field reaches 
saturation more slowly. 

3.4. Validity of our Approximation 

In Figure[3]the magnetic Prandtl number Pm is shown 
as a function of the density. Pm starts with an extremely 
high value of roughly 10^^ and then after a rather con- 
stant phase decreases rapidly. The magnetic Prandtl 
number is defined in equation (fTO|) with rj Tyohm + '7AD • 
For low densities the Ohmic resistivity dominates, which 
is independent of the magnetic field strength. With in- 
creasing density the magnetic field increases due to the 
dynamo amplification and with 77AD oc the ambipo- 
lar diffusion rate becomes dominant. In this regime the 



2 















ocr 






















oc £-'V4^ 










■s. 





















length scale £ 



Figure 5. Schematic model for the time evolution of the spectrum 
of the magnetic field strength in the inertial range of turbulence. 
For simplicity we use a fixed frame of reference, where the viscous 
and the Jeans scale stay constant. The different colors and line 
types represent the spectrum at different times. The green dashed 
line shows the spectrum at saturation on the viscous scale. To 
larger scales the spectrum decreases according to the Kazantsev 
slope {B oc £~^/4). At a time after saturation on the viscous scale 
the spectrum is indicated by the dotted red line. Finally, the solid 
orange line shows the spectrum, when the field is saturated on the 
Jeans scale. 

magnetic Prandtl number decreases proportional to B~^. 
As the magnetic field increases exponentially during the 
small-scale dynamo amplification in the beginning of the 
collapse, the magnetic Prandtl number decreases rapidly. 
The approximation of large magnetic Prandtl numbers 
(jSchober et al.l[20T^ is accurate during most of the dy- 
namo growth. At the end of the dynamo phase, however, 
Pm reaches unity and decreases even fur ther and our 
approximations eventually break down. Sc hober et al.l 
(j2012) show that for decreasing Pm the growth rate de- 
creases. However, they make no prediction for the regime 
Pm « 1. But numerical simulations sh ow that the dy- 
namo operates also in this regime (e.g. iFederrath et al.l 
l2011a[) . For Pm <^ 1 there is again analytic al evidence 
for s mall-scale dynamo action (e.g. Schckoc hihin et all 
l2007f ). We note that this treatment concerns the viscous 
scale only. 

4. MAGNETIC FIELD AMPLIFICATION ON LARGER 
SCALES 

4.1. Model for the Transport of Magnetic Energy to 
Larger Scales 

After the magnetic field saturates on the viscous scale 
the peak of the magnetic energy spectrum moves to larger 
scales. In this section we present a model for the time 
evolution of the magnetic energy spectrum. The situ- 
ation is illustrated schematically in Figure [5l Here we 
indicate three different curves, which represent different 
times. The dashed green line is the spectrum at the time 
of saturation on the viscous scale, the dotted orange line 
shows a later time and the solid red line represents an 
even later point in time at which the magnetic field has 
saturated on the Jeans scale. 

During saturation, the coherence length of the magnetic 
field shifts towards larger scales, a well-known behavior 
for the small-scale dynamo (|Schekochihin ct al.. ,2002bl : 
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Figure 6. Different cliaracteristic scales as a function of the den- 
sity. The dashed green Une indicates the viscous scale, the dotted 
red line the scale corresponding to the peak of the magnetic energy 
spectrum, and the solid orange line the Jeans scale. We show the 
results for Kolmogorov turbulence in the upper plot and the results 
for Burgers turbulence in the lower plot. 



iBrandenburg &: SubramaniarJ I2005D , r ecently shown t o 
be true also in a collapsing system (|Sur et al.l I2012D . 
Analytical arguments suggest this occurs on the eddy- 
timescale of the current peak scale £p 



IS 



(21) 



where we used i^p — vj{£p/£j)^. Considering that the 
peak scale moves from the viscous scale (-v{tu) towards 
larger scales, we find for the time-dependency of the peak 
scale 



^p(<)=4(t.)+(^(t-t.) 



(22) 



where is the point in time, when saturation occurs on 
the viscous scale. 

The slope of the curves proportional to is known as 

the Kazantsev slope in real spaccQ, which can be derived 
from the Fourier-tr ansformed Kazantsev equation (jl4p 
(|Brandenburg fc Subramanian ,2005,1 . This characteris- 

® In many references the magnetic energy spectrum is given as a 
function of the wave number k, defined for example as B'^ /(8np) = 
1/2 f M{k)dk. In this case the Kazantsev slope is M{k) oc fc''/". 
From this we find B" oc fc^/a and B oc £"5/4^ 



tic slope is also observ ed in simulations ()Federrath et al.l 
[mTilX^reFalllmil) . The curve that connects the peak 
maxima at different times (red-colored curve) is a relic 
of the turbulence spectrum and thus is proportional to 

At each time step we calculate the peak magnetic field 
strength by solving the stationary case of equation (I^Hl) . 
However, we find that the magnetic field strength ex- 
ceeds the equipartition field strength on scales larger 
than the viscous scale. The reason for this is that the am- 
bipolar dissipation rate, which is proportional to B'^/£'^, 
decreases rapidly in this regime. Thus, it cannot bal- 
ance the growth rate any longer and we need to set 
the equipartition field strength as an upper limit i?£.inax- 
With Bf^^^/{8TT) = l/2pv{tY we find the maximum 
magnetic field strength Bi^^i^^ — y/ATFpv{i). 

Taking the typical turbulent velocity on the scale of the 
turbulence £ to be related to the sound-speed by v{£) — 
{£/IjYcs ~ {-fkT/mf/^ii/ljY, we find that 



(23) 



Using the Kazantsev slope, we can extrapolate the mag- 
netic field strength onto the current Jeans length. By 
this we are able to determine the time evolution of the 
magnetic field on the Jeans scale. 

For this process to be relevant during collapse, the eddy- 
timescale needs to be smaller than the collapse timescale. 
Thus, the small-scale dynamo is unlikely to produce mag- 
netic fields on scales larger than the Jeans scale. Figure 
[B] shows the viscous, the peak, and the Jeans scale as 
a function of density. During the small-scale dynamo 
growth the spectrum of the magnetic energy peaks at 
the viscous scale. After saturation on the viscous scale 
the peak moves to larger scales according to equation 
(l22l) until it reaches the Jeans scale. 



4.2. Resulting Jeans-Scale Magnetic Field 

As described in the last section, we determine the mag- 
netic field on the Jeans scale by extrapolation from the 
peak scale. The result of the large-scale magnetic field 
is shown in Figure [7] together with the field on the cur- 
rent peak scale and the one on the viscous scale. One 
can see that the magnetic energy is shifted rapidly onto 
larger scales. For Kolmogorov turbulence the field on the 
Jeans scale saturates at a density of roughly 3 cm~^ and 
for Burgers at a density of roughly 4 cm~'^. At the end 
of dynamo growth on the Jeans scale we have a magnetic 
field strength of about 10^^ G throughout the entire iner- 
tial range of the turbulence, i.e. within the Jeans volume. 

After the rapid initial dynamo amplification the 
only way to amplify the magnetic field on the Jeans 
scale further is gravitational compression, which leads 
to S cx n^^^. However, the field has already reached 
equipartition with the kinetic energy at the end of the 
dynamo amplification and, thus, increases only with n^/^ 
(see equation [23| . The growth rate of the magnetic field 
on the Jeans scale Fj is then 



r.i 



1 dn 

n dt 



(24) 



In Figure |8] we compare the growth rate Fj to the am- 
bipolar and Ohmic diffusion rates on the Jeans scale. 
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Figure 7. The magnetic field strength as a function of the number 
density on different scales. The dashed green line corresponds to 
the field evolution on the viscous scale, the dotted red line to the 
peak scale and the solid orange line to the Jeans scale. We show 
the results for Kolmogorov turbulence in the upper plot and the 
results for Burgers turbulence in the lower plot. 




Figure 8. The growth rate on the Jeans scale Fj after the dynamo 
amplification compared to the diffusion rates as a function of the 
number density. Fohm.j a-nd Tad, J arc the Ohmic and ambipolar 
diffusion rate, respectively. 



Tad, J and Fohm.j- As Fj is always larger than the dif- 
fusion rates in the shown density range, the magnetic 
energy on the Jeans scale is not dissipated again during 
the collapse. At a density of 10^^ cm"^, we determine 
with B oc n^/^ a magnetic field strength of 0.4 G. 



4.3. Implications for Numerical Simulations 

Our calculations show that, due to the rather small vis- 
cosity and resistivity in primordial gas, the hydrodynam- 
ical Reynolds number, the magnetic Reynolds number 
and the magnetic Prandtl numbers have very high val- 
ues as long as the magnetic field is not saturated. Such 
Reynolds numbers are well above what can be reached in 
numerical simulations, implying that the physical growth 
rate of the magnetic field largely exceeds the growth rate 
obtained in numerical simulations. Particularly impor- 
tant here is the fact that the typically unresolved viscous 
scales are highly relevant for magnetic field amplification 
even on larger scales. In this sense, numerical simula- 
tions can only show the presence of a dynamo, but will 
typically underestimate the magnetic field amplification 
rate. This behavi or has also been dem onstrated in pio- 
neering studies bv lSur et all (|201Cl[ ) and iFederrath et all 
(IMIl . 

On the other hand, our results show that magnetic fields 
quickly saturate once turbulence forms, and the limit- 
ing timescale may thus be the timescale on which turbu- 
lence is generated. This is again an issue which can be 
addressed with numeric al simulatio n s, and indeed, sim- 
ulations for instance by iTurk et all (|2012D convincingly 
demonstrated the release of turbulence from the gravita- 
tional energy during primordial collapse. Overall, such 
simulations are thus relevant to explore the origin and 
generation of turbulence, while the strength of the mag- 
netic field should rather be estimated based on the phys- 
ical growth rates. As a net effect, we therefore expect 
that the magnetic energy is always close to saturation 
once turbulence is generated in the halo. 

5. SUMMARY 

We computed the evolution of the magnetic field and 
its saturation level in typical primordial halos based on 
the Kazantsev theory of the turbulent dynamo in com- 
bination with a detailed description of the physical and 
chemical processes in zero-metallicity gas. The model is 
in principle applicable only to magnetic field fiuctuations 
on very small scales. However, when interested in the in- 
fluence of the field on the overall dynamical evolution of 
the halo gas, it is most important to understand how sat- 
uration occurs on larger scales. To address this problem, 
we also considered the transport of magnetic energy from 
the viscous scale to the Jeans scale. 
Starting with a weak magnetic seed field of 10"^° G, as 
can be produced by the Biermann battery, we follow the 
evolution of magnetic field fluctuations on the viscous 
scale and found that they are amplified very rapidly on 
timescales much shorter than the free-fall time. As a con- 
sequence, the field saturates almost immediately after the 
onset of gravitational collapse in the halo. By extrapo- 
lating the small-scale magnetic field to larger scales and 
assuming the peak of the magnetic spectrum shifts on 
the local eddy timescale, we were able to follow the evo- 
lution of the magnetic field strength throughout the full 
inertial range within the Jeans volume. For typical halo 
parameters, the dynamo growth of the magnetic energy 
saturates at a density of roughly 3 cm""^ for Kolmogorov 
turbulence and 4 cm"'^ for Burgers turbulence. At this 
point in time the field has a strength of about 10"^ G. 
We point out, however, that the field continues to grow 
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in the collapsing gas due to gravitational compression. 
Our results show that the magnetic energy on small 
scales, and more importantly also on dynamically impor- 
tant large scales, can grow to very high values. In order 
to understand the influence of this strong field on the evo- 
lution of the halo gas, it is important to know whether 
the small-scale magnetic field can be transformed into 
a coherent large-scale field. One way to produce more 
coherent magne tic structur e s is b y forming disks, which 
is suggested bv iLatif et al.l (|2011| ). Moreover, the satu- 
ration behavior of the small-scale dynamo should be ex- 
plored further in the regime Pm < 1, as we have shown 
that the magnetic Prandtl number is in this regime for 
high densities. 

If indeed the processes discussed here can produce dy- 
namically significant fields on large scales, then the mag- 
netic field will infiuence the star formation process in 
high-redshift halos. For example, since recent high- 
resolution simulations indicate that the accretion disks 
around the very first stars were stro ngly susceptible 
to fragmentation (iTurk et all [2 009: St acy et all boid: 
IClark et al.llMl IGreif et al.ll20lli: iSmith et al.ll201llH t 
is expected that most primordial stars formed as mem- 
bers of binary or higher-order multiple systems with a 
wide range of masses rather than being isolated, high- 
mass stars. From studies of low-mass star formation at 
present day, however, we know that magnetic fields close 
to the equipartition value can effectively redistribute an- 
gular momentum via a process called magnetic braking 
(jMachida et al.ll2008bt iMouschovias fc Paleologou.19791) 
and can t hereby reduce the fragmentation probability in 
the disk (iHennebelle fc Ciardill200l iPeters et al.l 120111: 
iHennebelle et al.ll2011l: iSeifried et al.ll201lD . The "correct 
treatment of magnetic fields in calculations of primordial 
star formation therefore seems critical to better under- 
stand the mass function and multiplicity of metal-free 
stars. 
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